Implementation of logical gates on infinite dimensional quantum oscillators 
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Abstract — In this paper we study the error in the approx- 
imate simultaneous controllability of the bilinear Schrodinger 
equation. We provide estimates based on a tracking algorithm 
for general bilinear quantum systems and on the study of the 
finite dimensional Galerkin approximations for a particular 
class of quantum systems, weakly-coupled systems. We then 
present two physical examples: the perturbed quantum har- 
monic oscillator and the infinite potential well. 

I. Introduction 

A. Logical gates 

Quantum computation relies on the idea to store an infor- 
mation in the state of quantum system. This state is described 
by the wave function, that is, a point in the Hilbert sphere 
of C), where f2 is a Riemannian manifold. 

When submitted to an excitation by an external field {e.g. 
a laser), the time evolution of the wave function is governed 
by the bilinear Schrodinger equation 

dih 1 

^ ~-A4, + V{x)ij{x,t) + u{t)W{x)i:{x,t), (1) 

where V,W : H are real functions describing respec- 
tively the physical properties of the uncontrolled system and 
the external field, and u : R — > R is a real function of the 
time representing the intensity of the latter. 

When the manifold is compact, the linear operator 
i{A/2 — V) admits a set of eigenstates ((/>„)„gN-A logi- 
cal gate, or quantum gate, is a unitary transformation in 
L^(ri, C) for which some finite dimensional space of the 
form span{0i, 02, •■•,</>«} is stable. To build a given logical 
gate T from the system ([T]i, one has to find a control law 
u such that the propagator at a certain time T of ([T]) 
satisfies Tl^{(j)j) — T(j)j for every j = 1, . . . ,n. 

The main difficulty with this problem is that the space 
L^(ri, C) has infinite dimension. For the sake of simplicity, 
one often only considers the case where il is a finite union 
of points (or, equivalently, L^(n, C) is finite dimensional). 
Nevertheless, most of the usual quantum systems evolves on 
non trivial manifolds il. This papers deals with the effective 
implementation of some simple logical gates on models of 
quantum oscillators on 1-dimensional manifolds. 

B. Quantum control 

The problem of driving the solutions of ([TJ to a given 
target has been intensively studied in the past decades, both 



in the finite and infinite dimensional case. Many advances 
have been done in the infinite dimensional case, when there 
is only one source and one target. The interested reader 
may refer, for instance and among many other references, 
to [5], [7] for the theoretical viewpoint and to [8] for 
numerical aspects. In particular, it was proved in [3] that 
exact controllability is impossible in general. This does not 
prevent to study approximate controllability of Q, that is to 
replace the condition T1^{(f)j) — T(f)j by ||Ty((/)j) — T0j|| < 
e for every j = 1, . . . , n. To the best of our knowledge, there 
are only very few results of simultaneous controllability in 
the infinite dimensional case and the only available effective 
control techniques have been described in [9] and [11]. 

Recently, we noticed in [10] that a certain class of bi- 
linear systems are precisely approached by their Galerkin 
approximations. Two important examples of these so-called 
weakly-coupled systems are the quantum harmonic oscillator 
and the infinite potential well. The structure of weakly- 
coupled systems permits precise numerical simulations for 
the construction of quantum gates. 

C. Content of the paper 

The theoretical background is recalled in Section |ll] Be- 
sides a quick survey on simultaneous control techniques for 
Equation ([TJ, we give precise definitions and approximation 
results for weakly-coupled systems. In Section III we apply 



these results to a perturbation of the quantum harmonic 
oscillator and we provide estimates for the error in the 
controllability in a suitable finite dimensional approximation. 



Similarly, in Section IV we study the infinite potential well. 



II. General theoretical results 

A. Framework and notations 

We reformulate the problem ([T]) in a more abstract frame- 
work. In a separable Hilbert space H endowed with norm 
II • II and Hilbert product (•,•), we consider the evolution 
problem 



~dt 



= {A + uB)il) 



(2) 



where {A, B) satisfies Assumption [T| 

Assumption 1: {A, B) is a pair of linear operators such 
that 

1) for every u in R, A + uB is essentially skew-adjoint; 



2) A is skew-adjoint and has purely discrete spectrum 
(— iAfc)fegN, associated to the Hilbert basis {(f>k)keN 
of eigenvectors of A. 
From Assumption T[T one deduces that, for every piece- 
wise constant u, u : t I]^ MjX(t,,t,+i)(i), with = 
< ^1 < ■ • ■ < tN+i and uq, . . . , un in R, the solution 
1 1-^ T^^/jo of Q has the form 



o e 



(tj_l-tj_2)(^ + «3-2S) 



o e 



to(A+«o-B) 



^0 



for t e [tj_i, ij). 

Remark 1: With extra regularity hypotheses, it is possible 
to define the propagator of (|2]i for a larger class of controls. 
For instance, when B is bounded, for every t, Tj : u i— ^ TJ' 
admits a unique continuous extension to ^[^^(R, R). 

B. Control results 

Definition 1: Let (A, _B) satisfy Assumption [T| A subset 
S' of couples two levels j, fc in N, if there exists a finite 
sequence ((sj, Sj), . . . , (s^, Sj)) in S such that 

(0 '^1 = j and S2 = fc; 

(m) = s{^^ for every 1 < j < 9 — 1; 

M (<^,.,B(/.,,)7^0forl<j<<z. 

The subset is called a connectedness chain for (A, i?) if 
S couples every pair of levels in N. A connectedness chain 
is said to be non-resonant if for every (si, S2) in S, |Asi — 
AsJ 7^ I Ail -At J for every (^1,^2) in N2\{(si, S2), (s2,si)} 
such that {(j)t^,B(t)tJ ^ 0. 

Definition 2: Let S) satisfy Assumption [T| The sys- 
tem {A, B) is approximately simultaneously controllable if 
for every T e U{H) (unitary operators acting on H), 
. . . ,iljn G H, and e > 0, there exists a piecewise constant 
function Ug : [0, T^] R such that 

||TV;,-T^:V,II <£• 

for every j = 1, . . . ,n. 

The following sufficient condition for approximate simul- 
taneous controllability has been given in [9]. 

Proposition 1: Let [A, B) satisfy Assumption[r|and admit 
a non-resonant chain of connectedness. Then (A, B) is 
approximately simultaneously controllable. 

C. Weakly-coupled systems 

Definition 3: Let fc be a positive number and let {A, B) 
satisfy Assumption 1 1 Then {A, B) is k weakly-coupled 
if for every u e R, D{\A + uB\''/^) = D{\A\''/^) and 
there exists a constant C such that, for every ip in DdAI*^), 

|3?(|A|V,S/^>I <C|(|A|V»|. 
Definition 4: Let G N. We define the projection ttn ■ 
k<Ni'f'ki''P)4>k- The Galerkin approximation of 
(|2]i of order is the system in H 



where ^^^^ = ttnAttn and B^^^ = tttvBttjv are the 
compressions of A and 5 (respectively). 



We denote by X^j^s^ (t, s) the propagator of ( |S]jv[ ) associ- 
ated with a piecewise constant functions u. 

Proposition 2: Let k and s be non-negative numbers with 
< s < fc. Let {A, B) satisfy Assumption [T] and be k 
weakly-coupled Assume that there exist d>Q, Q<r<k 
such that < d\\i})\\^i2 for every V in D{\A\'/'^). Then 

for every e > 0, -fC > 0, n e N, and ('0j)i<j<n in 
D{\A\^'/'^Y'^ there exists G N such that for every piecewise 
constant function u 

MlkK =^ ||Tr(V^j)-^w(i,0Vjv^,|U/2 <£, 

for every t > and j = 1, . . . ,n. 

III. The perturbed quantum harmonic oscillator 

A. Physical model 

The quantum harmonic oscillator is one of the most 
studied quantum system. Schrodinger equation reads 



1 

'2 9a?2 



1 2 

-X 

2 



u{t)x j tl){x, t) 
R. With the notations of (|2]i, A = -i 



(3) 




where x £Vl - 
x'^)/2 and B — ix. 

An Hilbert basis of H made of eigenvectors of A is 
given by the sequence of the Hermite functions (0„)„gN, 
associated with the sequence (— jA„)„gN of eigenvalues 
where A„ = n — 1/2 for every n in N. In the basis ((^n)„gN, 
B admits a tri-diagonal structure: 

if j = fc - 1 
if j = fc + 1 
otherwise, 

A chain of connectedness for this system is given by S* = 
{{n,n +1) : n € N}. The chain S is resonant indeed 
|A„+i — A„| = 1 for every n in N. As a matter of fact, the 
system (|3]l is known to be non-controllable (see [14], [13]). 

We consider a perturbation of this system. Consider 
the inverse A~^ of the operator A. The family ((/)„)„gN 
is a family of eigenvectors for A^^ associated with the 
eigenvalues (— i/-^n)neN- For every 77 > we set = 
A + TjA^^. Since A and A^^ commute then {(f>„)neN is a 
family of eigenvectors for A,, associated with the eigenvalues 
(— iA5J)„gN where AJJ = A„ + rj/Xn- The set 5 is a non- 
resonant chain of connectedness for system (A^, B) for every 
7? > 0. Indeed A^^+i - A;^ = 1 - A^J^ and, clearly. 



X2 = a; 



V 

m+1 



AJ',, if and only if n — m. 



By Proposition [T] the system is approximately 

simultaneously controllable. Moreover by [9, Theorem 2.13] 
we have also an upper bound on the L^-norm of the control 
independent of the error. For instance we can steer approx- 
imately the first level (pi to the second 02 with a control 
law with L^-norm smaller than 57r/4. Another consequence 
is that a quantum gate for 0i,02, and ^3 is approximately 
reachable, that is for every e > 0, there exists > 
and a piecewise constant function such that \\T^^ {(f>j) — 
0CT(j)ll < £ where cr is the 3-cycle which exchanges 1, 2 and 
3. This can be achieved with < 7r/2(l + •\/2/2). 



B. Estimates 

In the following, we only consider control of _L^-norm 
less than K = 7r/2(l + \/2/2). The particular tri-diagonal 
structure of system {Ajj,B) is very useful for a priori 
estimates on the components of the propagator Indeed if 
< K, by [10, Remark 6], we have that 



(2if)"-2 /(2n-3)! 
(n-2)! V (n-2)! ' 



(4) 



for every n in N, n > 3 and j = 1, 2, 3. 

We use Q to find estimates on the size N of the Galerkin 
approximation whose existence is asserted by Proposition [2] 

First, let N > j and notice that 

|7r^Tr(</.,) = (AW+^i?W)7r^Tr(</),) 
Hence, by variation of constants, for every t > 0, 



+ / X("jv)(t,s)^ArB(Id-^Ar)T^(0,Xr)dT. (5) 

Therefore, since -'^(V) ^) unitary and for the tri-diagonal 
structure of the system we have, for j ~ 1,2,3, 

ik^Tr(0,)-X("^)(t,s)0,ii 



< 



< 



K\\7rNB(ld~7:N)r^{<j),)\\ 
K\bN,N+i\\{<l>N+urn<l>j))\ 




(2iV- 3) 



(TV -3)! ■ 

Using K = 7r/2(l + \/2/2) < 2.69, it is enough to 
consider a Galerkin approximation of size N = 420 to get 
||^jvTr (0,) - X^^^^it, < 10-4 for J = 1, 2, 3. 

C. Numerical simulations 

For simulations, we choose 77 = 1. To induce the transition 
between levels 1 and 2, the control law we use is a piecewise 
constant Aw periodic function, taking value 1 for < t < 
5. 10~^ and taking value for 5. 10^^ < t < 47r. We apply 
this control for 314 periods, that is during a time of 12567r . 

To induce the transition between levels 2 and 3, the control 
law we use is a piecewise constant 127r/5 periodic function, 
taking value 1 for < t < 5. lO^'^ and taking value for 
5. 10"'^ < t < 127r/5. We apply this control for 222 periods. 

The simulations are done on a Galerkin approximation 
of size 420, which garantees ||7rjvT"(0j) — X"^^(t, s)<j)j\\ < 
IQ-^ for j — 1,2,3. At final time T, the resulting propagator 
is such that |(T^0i,03)| > 0.998, |(T^(/.2, (/-i)! > 0.999 
and |(T^03, (j)2)\ > 0.999. The time evolution of the moduli 
of the first coordinates of T"(0j) for j = 1,2,3 is depicted 
in Figures 1, 2, and 3. 

All the computations were done using the free software 
NSP, see [12]. The source code for the simulation is available 
at [l].The total computation time is less than 4 minutes on 
a standard desktop computer. 




Fig. 1. Time evolution of the moduli of the first three coordinates of 
T"</)i in the case of the perturbed harmonic oscillator. First coordinate in 
blue, second coordinate in green, third coordinate in red. For the sake of 
readability, time scale is 1:100, total duration around 5500. 




Fig. 2. Time evolution of the moduli of the first three coordinates of Tf4'2 
in the case of the perturbed harmonic oscillator. First coordinate in blue, 
second coordinate in green, third coordinate in red. Time scale is 1:100. 



IV. Particle in a box 
A. Physical model 

We consider now the case of a particle confined in (0, tt). 
This model has been extensively studied by several authors 
in the last few years and was the first quantum system for 
which a positive controllability result has been obtained. 
Beauchard proved exact controllability in some dense subsets 
of using Coron's return method (see [4], [6] for a precise 
statement). Nersesyan obtained approximate controllability 
results using Lyapunov techniques. In the following, we 
extend these controllability results to simultaneous control- 
lability and provide some estimates of the L^-norm of the 
controls achieving simultaneous controllability. 

The Schrodinger equation writes 



dt 



1 d^tfj 

2 



— u{t)x'ijj{x, t) 



with boundary conditions tp{0,t) — ■ip{Tr,t) 
t e R. 



(6) 

for every 



so that 




Fig. 3. Time evolution of the moduli of the first three coordinates of T^ct>z 
in the case of the perturbed harmonic oscillator First coordinate in blue, 
second coordinate in green, third coordinate in red. Time scale is 1:100. 



In this case % — ((0, tt), C) endowed with the Hermi- 
tian product {^i,i^2} = Jq '^i{x)iJ2{x)dx. The operators 



A and B are defined by At/) 



'2 



for every ij: in 



D[A) = {H2 n Hi) ((0, tt), C), and B^j) = ixip. An Hilbert 
basis of H is {<j>k)ke'N with ipk ■ x t-^ sm{kx)/\^. For 
every k, A(j)k = —ik'^/2(j)k- 



For every j, k in N, 



bjk ■■= 




if j — fc odd 
otherwise. 



Despite numerous degenerate transitions, the system is 
approximately simultaneously controllable (see [9, Section 
7]). 

B. Estimates 

Using Proposition |2] to estimate the error done when 
replacing infinite dimensional system by its Galerkin approx- 
imation one finds, for = 97r/16 (see [10, Remark 4], 
with K = 97r/16, d = tt, fc = 1, r = 1, ci{A, B) < tt + 2, 
e = 10-3), that if TV > 1.6 10^ then 

|K^TJ'0i-X("^)(t,O)0i|| < 10-^ 

This estimation is definitely too rough to allow easy 
numerical simulations: matrix B^^^ ^ has about 5 10^^ non- 
zeros entries, the numerical simulations at such scale are 
difficult without large computing facilities. We have to go 
more into details to obtain finer estimates. 

Assume that, for some iV in N and 77 > 0, the control 
M : [0, T] ^ R is such that, for every t in [0, T], 



|X("^)(i,0)^3-^3^(V)(<,0)|| <r;. 



We have 



7r3X5v)(t,s)-^(V)(i, 5)713 

= ^3^(iv) it, 0)X('V) (0, s) - (t, 0)X("^) (0, s)^3 
= 0)(7r3X(V)(0, s) - X("^)(0, 5)^3) 

+(7r3X(V) {t, 0) - [t, 0)^3)^("w) (0, s) 



||7r3X(V)(i,5) - 5)71311 < 2r,. (7) 

Projecting (|5]) on the first 3 components we have, for j = 
1,2,3 that 

IK3Tr(</',)-^3Xw(t,o).^,l! 



< 



< 



|7r3X("^)(t,s)^ArS(Id - 7rAr)T,"(0,)||u(s)ds 

/ ||X("^)(i,s)7r3B(Id-7rAr)T:(<^,)h(s)ds 
Jo 

+ I ||(^3^("^)(i,s)-^rJV)(i,s)^3)|i||S||u(s)ds 

\u{t)\d?j (|K3B(Id-^Ar)|l 
+2||S||sup||7r3X("^)(t,0)-X('V)(i,0)^3)ll).(8) 



By skew-adjointness, ||7r3_B(Id— 7rAr)|| — ||(Id— 7rjv)i37r3||. 
This last quantity tends to zero, and we are able to give 
estimates of the convergence rate. Indeed, 

2 



< 



\{U-ttn)B<Pi 



< 



< 



< 



E 

k>N 

k>N 
1 



2k 



(fc- 1)2(1 + fc)2 
1 



(k-iy 



Similarly, 



i(Id-7rAr 



< 



< 



V2 

(iV- 3)5 



(iV-4)5' 

The procedure to induce a given transformation, up to a 
given tolerance e > 0, on the space span{0i, 02: ^3} is the 
following: 

1) Use estimates given in [9] to give an a priori upper 
bound K on the L^-norm of the controls one will use. 

2) From K and e, find N such that ii:||7r3B(Id-7rAr)|| < 
e/2. 

3) In the finite dimensional space span{0i, . . . , ^jy}. 
consider the bilinear system i = (A^^^ + uB'^^^)x 
and find a control u achieving the desired transition 
up to e/{2K) and such that < K. This can 
be done using standard averaging procedures (see for 
instance [15]). 

4) Use (|8]l to get an upper bound of the distance of the 
trajectories of (Sat) and the actual infinite dimensional 
system. 

C. Numerical simulations 

We illustrate the above procedure on an example. Fix 
e = 7 10-2. We would like to find u : [0, T] ^ R such 

that |(</)3,T^</>i)| > 1 - £, |(</.i,T^(/.2)| > 1 - £ and 




20 40 60 80 100 120 




Fig. 4. Time evolution of the moduli of the first three coordinates of T^fpi 
in the case of the potential well. First coordinate in blue, second coordinate 
in green, third coordinate in red. 



Fig. 5. Time evolution of the moduli of the first three coordinates of T"(/)2 
in the case of the potential well. First coordinate in blue, second coordinate 
in green, third coordinate in red. 



\{(j)2, '^t4'3)\ > 1 — £ at final time T. For this example, we 
are not interested in the respective phases but the method can 



easily be generalized to address this point (see Section IV-D 
below). 

From [9], the transition can be achieved with controls of 
L^-norm smaller than 57r/4(9/8 + 25/24). Using controls 
with better efficiencies (as described in [11]), we can use 
controls with L^-norm smaller than 2(9/8 + 25/4) = 13/3. 

Using the above estimates, one sees that if 20, then 



i^|j7r3B(ld-7rAr)|| < 



13 



< 6 10" 



3 (TV -4)5/2 

Last, we define u by u{t) = cos(3f)/20 for < t < 
72 and u{t) = cos(5i)/20 for 72 < t < T = 138. We 
check that J^^ \u{t)\dt < 13/3. One checks numerically that 

|k3^("20.(i,0) -X." (i,0)^3li < 1-3 10^' for t < 138. 



From (|7]), we get, for every t,s < T 

||^3X(V)(i,s) -X5v)(t,s)^3|| < 2.6 10-3. 
From ([8]l, we have, for j = 1,2,3, 

||7r3TJ'(0,)-^3^('V)(i, 0)0,11 <y(6 10-3 + 8.2 IQ-^) 

< 6.1 IQ-^ 



Conclusion follows from the numerical computations 

|(03,^("2o)(I^,O)</.i)| « 0.99924 

|((/.i,X("2o)(T,O)02)| ~ 0.99943 

I (02, ^("20) (T, 0)03) I « 0.99949. 

The actual precision is likely much better than 6.1 10^^ 
which is known for sure. However, our estimates do not allow 
a better conclusion. 

The evolutions with respect to the time of the moduli of the 



first coordinates of Xji,m(j)k for k = 1,2,3 are represented 



in Figures 4, 5, and 6. 

All the computations were done using the free software 
NSP, see [12]. The source code for the simulation is available 
at [2]. The total computation time is less than 4 minutes on 
a standard desktop computer. 
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Fig. 6. Time evolution of the moduli of the first three coordinates of Tf4>3 
in the case of the potential well. First coordinate in blue, second coordinate 
in green, third coordinate in red. 



D. Possible improvements 

If one is interested not only in the modulus but also in the 
respective phases of the final points, it is enough to replace 
the functions 1 1-^ cos(3t)/20 and 1 1-^ cos(5f)/20 above by 
t H> cos(3t + 9i)/20 and 1 1-> cos(5i + 6'2)/20 respectively, 
where 61 and 6*2 are suitable phases. 

In order to get better precision in the approximation (i.e. 
a smaller e), it is enough to replace the functions t 
cos(3i)/20 and t M- cos(5t)/20 above by the functions 
t i-> cos{3t)/L and t 1— > cos{5t)/L with L large enough. 
The price to pay for a better precision is an increase in the 
time needed for the transfer. 

V. Conclusion and future works 

We have shown how it was possible to implement a 
quantum gate on two types of infinite dimensional quantum 
oscillators. Our method provides rigorous estimates and 
permits numerical simulations that can be run on standard 
desktop computers. 

A limitation of our models is that the Schrodinger equation 
neglects decoherence. This approximation may be justified 



for time small with respect to the relaxation time of the quan- 
tum system. Future works may focus on the optimization of 
the time of implementation. 
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